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ABSTRACT 

Reynolds stress closure models based on the recursion renormalization group theory are devel- 
oped for the prediction of turbulent separated flows. The proposed model uses a finite wavenum- 
ber truncation scheme to account for the spectral distribution of energy. In particular, the model 
incorporates effects of both local and nonlocal interactions. The nonlocal interactions are shown 
to yield a contribution identical to that from the e -RNG, while the local interactions introduce 
higher order dispersive effects. A formal analysis of the model is presented and its ability to accu- 
rately predict separated flows is analyzed from a combined theoretical and computational stand 
point. Turbulent flow past a backward facing step is chosen as a test case and the results obtained 
based on detailed computations demonstrate that the proposed recursion -RNG model with finite 
cut-off wavenumber can yield very good predictions for the backstep problem. 
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1. INTRODUCTION 


Turbulent flows of scientific and engineering importance are characterized by a broad spectrum 
of length and time scales. While the physical aspects of turbulent flows are best described by the 
equations of motion, limitations in computer capacity and speed preclude their direct solution for 
complex flows of relevance to technical applications. The current practice for high Reynolds num- 
ber flows of practical interest therefore involves some type of modeling for Reynolds stresses. The 
commonly used turbulence models are based on the calculation of one-point first and second mo- 
ments such as the mean velocity, mean pressure and turbulent kinetic energy. Among these, the 
two-equation turbulence models that involve the use of transport equations for the turbulent field 
parameters that involve the length and the time scales are probably the most widely used. They in- 
volve the simplest level of Reynolds stress closure that do not depend specifically on the flow ge- 
ometry. (For an excellent review of recent trends in analytical methods for Reynolds stress closure, 
the reader is referred to Speziale. 1 ) 

In its standard form the two-equation Reynolds stress turbulence models involve the turbulence 
kinetic energy and dissipation based on a Boussinesq type approximation of the form 


^ = -i*V v r 


'dU i dUj 


wherein U is the mean velocity based on Reynolds average, K is the turbulence kinetic energy, and 
Vyds the eddy viscosity which is isotropic. Such a representation of turbulence is often not effective 
from both theoretical as well as phenomenological point of view and the shortcomings associated 
with it are fully discussed by Speziale. 1 To over come some of these, models that are nonlinear 
(i.e., quadratic) in the mean strain rate were proposed in the form of a constitutive relation. 3,4 Spe- 
ziale 3 employed tensor and dimensional analysis, together with invariance constraints, to derive his 
model. Yoshizawa’s model 4 was obtained by appealing to a two-scale direct interaction approxi- 
mation. The application of these models depend on the empirical evaluation of the model constants. 
Apart from the specific values of the constants in these models, these two models have quite similar 
structure and both were able to predict the anisotropy in the Reynolds stresses in a noncircular duct 
problem. 

The present study addresses the need for a more effective approach in the development of two- 
equation turbulence models, and in this context the renormalization group ( RNG ) theory based 
models are examined for further development. While the application of renormalization group the- 
ory to turbulence has attracted much attention, 5 ' 18 it is important to realize that these calculations 
fall into two distinct categories: (a) e-RNG 5 ' 9 , pioneered by Forster et. al. 14 , and (b) recursion- 
/?A/G, 10 ' 13 pioneered by Rose. 15 These techniques have been critiqued 6,16 ' 18 and compared. 12 
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Here, we wish to point out that in the € -RNG, a small parameter e is introduced into the exponent 
of the forcing correlation function (with the forcing function being introduced into the momentum 
equation). The theory is then developed for e « 1, and all constants generated are evaluated in 
this limit e « 1. However, at the same time, all exponents that are e -dependent are evaluated at 
e = 4. 17 In fact, e = 4 is required in the e -RNG to recover the Kolmogorov energy spectrum in 
the inertial range. 17 This not only plays some havoc with the evaluation of constants, it also leads 
to another problem: /?NG-induced interactions that can be shown to be irrelevant in the limit 
€ « 1 cannot now be shown to be irrelevant in the limit € — > 4. Yet in e -RNG theories (which 
require e — » 4) these higher order nonlinearities are assumed to be unimportant. 6,12,16 Moreover, 
e - RNG theory can only take into account non-local interactions. 6,10,16 However, the procedure 
is quite amenable and Rubinstein and Barton 7 have derived a Reynolds stress model using e -RNG 
methods that is qualitatively similar to that of Speziale 3 and Yoshizawa. 4 

On the other hand, recursion-/?NG does not rely on an e -expansion, and treats explicitly the 
cubic nonlinearities induced into the renormalized momentum equation. Moreover, recursion-/? A/G 
can handle both local and non-local interactions. Effects such as the cusp behavior in the transport 
coefficients as k — » k c are recovered in these theories (here, k c is the cutoff wavenumber separating 
the large scale from the local resolvable scales). 13,19 ' 24 These effects are the consequences of local 
interactions and the cubic nonlinearities introduced by the RNG procedure. However, one of the 
major difficulties to the application of recursion-/? A/G for turbulent flows governed by the Navier- 
Stokes equation was that the eddy viscosity was determined as a fixed point of a very complicated 
integro-difference equation. This drawback has now been removed by extending the theory to han- 
dle the iterative removal of infinitesimal wavenumber bands. 25 Now, as in e -RNG, the eddy vis- 
cosity is readily determined from the solution of a relatively simple differential equation. Unlike 
e -RNG, though, the transport coefficients are determined over the whole resolvable scales and not 
just in the wavenumber limit k — » 0. 

In the present work, the recursion RNG procedure is used to develope a Reynolds stress closure 
model in a formal manner. It is shown that the recursion-/? A/G based model introduces two addi- 
tional terms arising from the local interactions effects that are of higher order than those obtained 
using the e -RNG method by considering only the long wavelength, non-local interaction limit. The 
first of these two terms as well as those from the conventional turbulence models are shown to be 
a part of the integrity basis used in the representation of the anisotropic part of the Reynolds stress 
tensor. The second term which arises from pressure-strain coupling is quartic in strain rate and 
is of the same order as the remaining terms of the integrity basis. 

The model is then applied for turbulent flow past a backward-facing step which has played a 
central role in benchmarking the performance of turbulence models for separated flows. During the 
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past decade — beginning with the 1980/81 Stanford conference on complex turbulent flows 29 — 
various two-equation turbulence models have been tested and compared with the experimental data 
of Kim, Kline and Johnston 30 and Eaton and Johnston 31 for the backstep problem. Initial results 29 
indicated that the standard K-z model, with wall functions, underpredicted the reattachment point 
by a substantial amount on the order of 20-25%. Several independent studies have been subse- 
quently published using alternative forms of the K-z model wherein a variety of conflicting results 
have been reported. Considering the need to accurately predict separated turbulent flows — which 
can have a wealth of important scientific and engineering applications — the proposed model is 
applied for the backstep problem. The computations based on a sufficiently resolved finite-volume 
algorithm show that the proposed model based on the recursive application of the renormalization 
group theory (developed independently without any ad hoc empiricism) can yield a prediction for 
the reattachment point that is within a few percent of the experimental result. The physical impli- 
cations that these results have will be also discussed in detail in the sections to follow. 


2. FORMULATION OF THE PHYSICAL PROBLEM 


The turbulent motion of viscous, incompressible fluids are governed by the Navier-Stokes 
equations which may be analyzed using classical single point closure based on Reynolds decom- 
position of all physical variables. The resulting averaged equations of motion are of the form 


W, 4 . n dU ‘ = -’ >p , * U ‘ _ 

a dx a 5^ V °dx a dx a dx a 


Wa 

dx a 


= 0 


( 1 ) 

( 2 ) 


where £/, is the mean velocity, P is the mean pressure, v 0 is the kinematic viscosity of the fluid and 
%ij is the Reynolds stress tensor. While the physical aspects of turbulent flows are best described 
by the above equations, limitations in computer capacity and speed preclude their direct solution 
for complex flows of engineering importance. The current practice for high Reynolds number 
flows of engineering interest therefore involves some type of modeling for Reynolds stresses. The 
commonly used turbulence models are based on the calculation of one-point first and second mo- 
ments such as the mean velocity, mean pressure and turbulent kinetic energy. 

In the present work we consider the development of Reynolds stress model by the recursion 
RNG formulation and its application. In this context it is convenient to express the equations of mo- 
tion using the Fourier representation, 


( 3 } + V 2 ) «, (k> 0 = Wfajj (*) jd 3 P «a (p. 0 “p (k - P, 0 


( 3 ) 


k a u a (k,t) = 0 ( 4 ) 

It is important to note here that no random forcing is introduced unlike in the e -RNG theories 
where it plays a critical role in introducing the small parameter € . The nonlinear coupling coeffi- 
cient 


= kpD ia (k)+k a D ifi (k ) 
with D a p being the projection operator defined by 

k ka 


( 5 ) 


( 6 ) 


A scale factor h is now introduced to partition the wavenumber space to N segments such that 
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( 7 ) 


^c~^N~ • 

kjij _ j = h ko , • • • > k j hk 0 

k 0 =O (k d ) 

where, k c is the cut-off wavenumber that separates the large scale from the small scale, and k^ is of 
the order of Kolmogorov dissipation wavenumber. The velocity field is then decomposed in the 
Fourier space such that u,- = uf 0 (k - k h ) + u*9(k h - k) (wherein 6 is the Heavyside unit step 
function, */,(/,= i;v) is the local cut-off wavenumber, the superscript < represents large resolvable 
scale quantities while > represents small unresolved scales with respect to the local cut-off wave- 
number). The evolution of these components can be directly obtained from (3) as: 

+ V °* 2) U>i (k ’ 0 = M,a P {k) r P {U<a (P ’ 0 “P (k " P ’ l) (8) 

+ 2 ul (p, t ) wj* (k - p, t) + (p, 0 (k - p, r)] 

+ v 0 * 2 ) uf (k, t) = M ia p (*) jd 3 p [ul (p, t) uj (k - p, r) 

+ 2 u> (p, t ) «p (k - p, t) + (p, t ) «p (k - p, r)] 

It should be noted that in e -RNG approach, one is forced into taking the large-scale infrared 
limit it -> 0. In essence, this forces a spectral gap between the resolvable part of the flow field and 
the small unresolved scales. If this spectral gap were somehow present initially, it would be quickly 
populated in just a few eddy turn over times. Thus, retaining only the distant interactions may not 
be appropriate. In fact, it has been shown 13 that the energy transfer function that corresponds to 
local interactions accounts for most of the energy flow out of the resolvable scales. It thus seems 
important to retain both local and nonlocal interactions in the modeling of the Reynolds stress and 
this can be readily achieved by recursion-/? A'G. In particular, it is apparent that the Reynolds stress, 
% = x^? + , has two components. The x^— part arises from the infrared limit of k — » 0 and 

is due to the u>- distant interaction limit while the - part arises from the 0 < k < k c spec- 
trum and is due to the uf - uf local interaction limit. Thus in the e - RNG model, the Reynolds 
stress X;j = r** and is obtained purely from the - wf interaction in the small unresolved scale 
momentum equation (8). We now consider the contribution to the Reynolds stress that arises from 
the local interaction: 

x>< =-jd 3 p[uf(k-p)uJ(p)+uf(k-p)uUp)] (10) 

Herein u < (k-p) corresponds to the Fourier-transformed velocity field in the resolvable large scales, 
Ik-pl < fcj, while u > (p) corresponds to the small scale field with Ipl > k j and k l is the wavenumber 
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which separates the resolvable from the small unresolved scales. 

Now there are certain constraints that X^ < must satisfy; in particular, 

(a) ^ 0 as the turbulent'kinetic energy K —> 0. 

(b) T ><: is Galilean invariant. 

Constraint (a) is clear: as the turbulent kinetic energy K — > 0, the small scale velocity field iv* — > 0, 
so that — > 0. Constraint (b) arises since it has been shown in detail that recursion-RNG, even 

with its cubic nonlinearity interaction, is a Galilean invariant theory. 25 Now, if one follows stan- 
dard recursion-/?NG procedures, 10 ' 13 the relevant part of the small scale velocity field that contrib- 
utes to xj^ is (on the removal of the first subgrid shell, c.f., equation (9)). 

(P) = — j d3 J «a(j) M p (j — P) + ••• (11) 

V oP 

where, . . . refers to terms that will not contribute to X^ . At the removal of the first small scale shell 
V 0 is the molecular viscosity. At subsequent steps in the removal of the subgrid scales, v Q will be 
replaced by the renormalized eddy viscosity. 

In the substitution of the subgrid velocity field from (11), into the Reynolds stress x^ given 
by (10), one must ensure that the constraints, (a) and (b) above arc satisfied. It can be shown that 
these constraints are not satisfied by a simple substitution of (11) into (10). Hence we introduce a 
factor (p - k) a (p - s ) a /k 2 c so as to recover these properties. This factor is appropriate since, 
from (11), k c = 0(e/K^ 2 ) — » °o as K — » 0. Hence, in recursion-/?NG, after the elimination of the 
first spectral band (i.e., after the first iteration) (10) becomes 


T ><: (k) = 

y v ' i 


J-w 


( Pa. kg) ( P a Sg) 

VqP 2 


Mjy 8 ip) Wy (5) «8 (p - s) u f (k-p) + i 


( 12 ) 


where p is in the subgrid region, while the arguments of all the velocity fields are in the resolvable 
scales. The second term in (12) is obtained by i <-» j interchange in the first term. At the next ap- 
plication of the above procedure (i.e., after the second iteration) 


tj* (k) l 2 = tJ < (k)| i -Jrf 3 /j[«<(k-p)M ; > (p) + uj (k — p) (p) ] (13) 

where, k < k 2 . For k 2 <k<k\, the small scales are obtained by a similar procedure by utilizing the 
equations of motion for small scales. 10,15 The resulting expression for the small scales may then 
be written as: 


«f <P> 


^iaP f ,3 • 

Vj (p)p 2 * J 


Ug (j)Mp (j-P) + ••• 


(14) 
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substituting from (12) and (14) in (13), 


l iT (k) ■ -p X [K^ 3 


C Pa k a ) (Pa 


y» = 1L 


V A (p)p 


2 A^/yg (P) My ( j) m| (p - S) uf (k - p) + I j 


(15) 


Proceeding similarly to the elimination of the N-th spectral band, one finds 

(P<x ~ k a) (Pa ~ s a) / \ < / \ < / \ < / / \ 

~ M jy 5 (P) M y (*) U S (P ~ s ) u i (k-p) + 1 


*r< k > =-pX[pv 


c A = 1 L 


V h (P)P 


(16) 


where, & c is the cut-off wavenumber that separates the large scales from the small scales. In (16) 
the summation is over all the spectral bands eliminated and V/,(p) is the renormalized eddy viscos- 
ity .10*13 now proceed to the differential limit of infinitesimal spectral shells, so that 1 ^ 


IrVP) ~*v(k c ) [*J 


4/3 


(17) 


for a Kolmogorov k inertial range energy spectrum. Since the quadratic coupling coefficient 


M jyd (p) = ^[ 5 ;yP» + Vy“ 2 ~?) 


(18) 


in the differential limit for infinitesimal spectral shells of the small scales, (16) becomes 


x>< (k) = 

y 2iv(k c )k l c 




■X X (Pa~ k a) (Pa~ s a) 

d 3 pd 2 s yti (Pp^ (s)up(p-s)u i ( k-p ) 


PqP jPy 

+ PpMp (s)Uj (p-s)uf (k-p) -2——^ — «p (s) Uy (p — s) uf (k-p) } +f«->y 

(19) 

The inverse transform of (19) yields the contribution to Reynolds stress due to local interaction in 
physical space and is shown below. 


*7? (*) = 


l 


(2lt) 3 V(* U 


10/3 


a /dt/ Y (x’)a 2 t/p(x’)> l 

+^ 2 (x x,) a?:[ a^ p dx' a dx' y J } 


( 20 ) 


In the above, the resolvable velocity field in the physical space is identified with the mean velocity 
Uj, and the nonlocal kernels gj(x - x') and # 2 ( x - x') are defined by: 


- 7 - 


r exp (ip ■ r) 3 , . r 

£i(r)=J 273 dp and 5 2 ( r ) = J 


f exp (ip • r) 3 


, 8/3 


d 3 p 


( 21 ) 


The structure of these nonlocal kernels is dependent on the cut-off wavenumber, k c . Now, 


gi(r) = ^ JV /3 sin (pr)dp = y k 4 / 3 jp' 1/3 sin (k c p’r) dp' 


2K 


4 k . 


( 22 ) 


where the integration limits are obtained from the partial average 15 over the nearest unresolved 
band. Similarly, 


g 2 (r) = ^k- c 2/3 jp'~ 5/3 sin (k c p'r) dp'. (23) 

l 

These nonlocal kernels reduce to local expressions for k c » 1 . Noting that sin (k^) /x — > k5 (x) 
as k c — » oo, 

g 1 (r)/4^ = 3.58^ /3 5(r) 

g 2 (r) /4* - 1.74*; 2/3 5 (r) . (24) 

Following Yakhot and Orszag, 5 this cutoff wavenumber k c can be related to the isotropic part of 
the turbulent kinetic energy K as: 


2/3 

Jf = j E (p)dp = lc* (25) 

k c k c 

on using the Kolmogorov intertial range spectrum for E(p), and where C K is the Kolmogorov con- 
stant. Thus, the cut-off wavenumber 


K = 


572 ^ 


3/2 


K 


Hence under this approximation xjjj 5 (x) in (20) reduces to a simple algebraic form 


(25) 




dUidUjBUp ; 


+ c 


K' 

R2~5 
£ J 


du i a (dUp d 2 u. 




dx„ l dx„dx n J + 1 




(27) 


where use has been made of Kraichnan’s 32 result 

v(k c )k 4 / 3 = 0.19C|e 1/3 


(28) 
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However, for most flows of interest, k c varies considerably throughout the flow domain. Thus, if 
one wishes to use the local approximation for the kernels, gi and g 2 as well as retain k c ~ 0(1) 
effects so as to obtain an algebraic form for (x) , then the coefficients C R1 and C R2 will be 
function of the flow quantities. In particular, C R i and C R2 may be expressed as: 

= C/jj (e, K, tj, £) and C R2 — C R2 ( £ > "H* 0 (29) 


where 


x fdu a bu^\ 

the strain rate, r| = (S ap S a p) 1/2 , with S a p ■ ^ + J and 

l (dU a BU$\ 

the rotation rate, £ ■ (^ap^ap) 1/2 > with W ap “ 2 ^9^ ” ~dx~ )' 


l P 

Rubinstein and Barton 7 have derived a Reynolds stress model using the e -RNG method (which 
corresponds to the infrared limit of k -> 0) of the following form: 



5 *V v r 


fdUi dUA K 3 

’5*1 J 




xl 


a ^ x aj 


pugv. W?Ez\\ c 

T2 dx ) d x a J t3 


fdU a dU a ' 

$ x i tej. 


(30) 


The constants C zl = 0.034, C x2 = 0. 104 and C T 3 = -0.014, and (. . .)* denotes the deviatoric part of 
the expression within the parenthesis. The first two terms correspond to the linear model and 
Vj. = C^K 2 /e is the isotropic eddy-viscosity, where £ is the turbulence dissipation and » 0.09 
based on empirical data from equilibrium boundary layer flows. The above model is quadratic in 
mean strain rate, includes the effect of convection and diffusion and is qualitatively similar to other 
second order models. 3,4 


Combining (27) - (30) a formal expression for Reynolds stress which includes both the local 
and nonlocal interactions may be expressed in the following form: 




(BU t BUA K 2 


<BU ; dU i 


V 




(BU t BU a BUjBU^ (BUa*U a >T 
x2 [dx a %Xj ) x3 ( Bx t ~5x~ J 


+ C, 


dU:dU:dU* 


K* 

R1 e 3 l dx a ^p^ 


. r K 
+ c /?2~5 

e 


(31) 


'du t a fBU} d 2 u y \ 

^x a ^Xj ^ dxy Bx a Bx^ j 


+ i <r*j 
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It is of some interest to recast the above expression for the Reynolds stress obtained from the 

recursion RNG model into an integrity basis representation which are commonly employed to rep- 

2 

resent the anisotropic part of the Reynolds stress tensor, b • = . - ^AT5. ) /2 K in the following 

form-invariant manner. ° 

b ij = (32) 

where, T,'- is the integrity basis for functions of symmetric and antisymmetric tensor and the co- 
efficients G are scalar functions of the irreducible invariants of the strain rate tensor, and the 
rotation rate tensor, Vf i y. 26 ' 28 

T ( i) _ o 

7'( 2 ) - f w v 

1 ij ia " a; " ia^aj 

T fP = C C . _ i C 5 

7 y °<a J oy ^ J aa u ij 

4 6> =• + 

7y S) = 

4 9> = W , »M'a P '5pA + S i „SafV 1 ,-|(S„DV^M' £0 )S i> 

r,)' 0) = W , to S 0 pS pT H' ye »'y-H' to »' op S^ 1E Vl'y (33) 

r\ n 

Recently, Gatski and Speziale have applied this integrity basis to determine an explicit alge- 
braic stress model for three dimensional turbulent flows based on a systematic derivation from a 
hierarchy of second-order closure models. In their applications, they have restricted themselves to 
a model 28 which involves just the tensors , and Tfp . It can readily be shown that the 

k = 0 part of the Reynolds stress, , specified in (30), involves the tensors , T-^ , T-p and 
T;P from the integrity basis given by (33). The X -part of the Reynolds stress is common to both 
e -RNG and recursion-RNG. However, unlike e -RNG, there is now a finite-/: spectral contribution, 
to the Reynolds stress in recursion-/?NG. It can be shown that the first term in (29) involves 
the tensors T-j ’ and T-- ’ of the integrity basis (33). The second term in (29), however, is not 
readily expressed in this second order basis because of the intrinsic fourth rank tensors involved. 
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However, as can be seen, this term is quartic in strain rate and arises from pressure-strain coupling, 

(7) (8") 

and could be thought of as having contributions similar to those involving the tensors T-j , T-j , 
and T'P . 

The above expression for Reynolds stress (31) are to be used along with the equations of mo- 
tion (1) by specifying turbulent kinetic energy and dissipation. In two-equation turbulence models, 
this closure is achieved through the development of transport equations for the turbulent kinetic 
energy and dissipation — quantities that are directly related to the length and time scales — of the 
following general form: 




T-z + 


dx : 



a/r 

ax- 


04) 


?l + u— = C -<P-C - + — 
dt U JdX: L ^K t2 K dX: 


( v 70 ae' 

i v ° + ^K 


% 


(35) 


where, v = Vq + v T is the total viscosity, IP = -Xy (dU/dxj) is the turbulence production, e is the 
scalar turbulent dissipation rate, and 2v 0 S l - / - ( duj/dx ; ) ( du/dxj ) is the turbulent strain rate 
correlation term. The quantities C t j, C&, a K , a e are dimensionless and taken to be 1.44, 1.92, 1.0 
and 1.3, respectively, consistent with the standard form of the two-equation K-z model (based on 
empirical data obtained from equilibrium boundary layer flows). 4 It should be noted here that these 
are calculated explicitly to be 1.42, 1.68, 0.719 and 0.719, respectively, in the e -PNG-based for- 
mulations. 5,6 In addition, it is customary to neglect the strain rate correlation term !Rjn the standard 
K-z model as well as in some PNG-based models 7 since 3£must be modeled in order to achieve 
closure. Special forms of approximations have been recently employed in the case of some of the 
e -PNG-based models. 8,9 To avoid ambiguity (and to avoid introducing ad hoc assumptions re- 
quired to model the contribution from the turbulent strain rate correlation term is neglected in 
the present work. 
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3. TURBULENT FLOW PAST A BACKWARD-FACING STEP — A CASE STUDY 

The problem to be considered is the fully-developed turbulent flow of an incompressible vis- 
cous fluid past a backward-facing step (a schematic is provided in figure 1). Calculations will be 
conducted for an expansion ratio (step height: outlet channel height) E of 1:3 and the Reynolds 
number Re = 132,000 based on the inlet centerline mean velocity and outlet channel height (which 
corresponds to that of Kim et a/. 30 and Eaton and Johnston 31 ). The mean turbulence equations (cf., 
§2) are solved subject to the following boundary conditions: 

(a) inlet profiles for U, K and e are specified five step heights upstream of the step comer ( U is 

™ 3 1 

taken from the experimental data ’ and the corresponding profiles for K and e are computed 
from the model formulated for channel flow), 

(b) The law of the wall is used at the upper and the lower walls, and 

(c) Conservative extrapolated outflow conditions are applied thirty step heights downstream of the 
step comer; these conditions involve the following: i) the V-component of the velocity for the 
cells at the outflow boundary are obtained by extrapolation; ii) the {/-component of the velocity 
is then computed by the application of a mass balance; and iii) the scalar quantities such as pres- 
sure, turbulent kinetic energy and turbulent dissipation are all obtained by extrapolation. It was 
found that a downstream channel length of about thirty step heights was needed to ensure that 
the local error for all the quantities was of the same order as the interior values. 

A finite volume method which relies on solving the discretized equations by a line relaxation 
method with the repeated application of the tridiagonal matrix solution algorithm is modified and 
applied for the present case to obtain the steady state solution. 33,34 The computed solution was as- 
sumed to have converged to its steady state when the root mean square of the average difference 
between successive iterations was less than 10' 4 for the mass source. 34 Approximately 2000 itera- 
tions were needed for the convergence of the standard K-e model; this corresponds to approximate- 
ly 20 minutes of CPU time in a partially vectorized mode on the Cray-YMP supercomputer using 
64-bit precision. The recursive RAC-based K-e model requires approximately 33% more CPU time 
due to the fact that the additional terms in the RAG-based K-e model have to be evaluated during 
each iteration. These correction terms are dispersive — an additional feature that slows conver- 
gence. 

The issue of resolution is crucial for the backstep problem and calculations indicate that a 
200x100 mesh yields a fully grid independent solution. 33 All of the computations conducted in this 
study were performed using this 200x100 nonuniform mesh. As indicated earlier, the inlet condi- 
tions were specified 5 step heights upstream of the step comer and the outlet boundary conditions 
were specified 30 step heights downstream of the step corner. It is crucial that a sufficient distance 
downstream of the reattachment point be allowed before the outflow conditions are imposed. Many 
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earlier computations of the backstep problem were in significant error due to the imposition of ful- 
ly-developed outflow conditions too close to the reattachment point. Furthermore, it is crucial that 
a fine mesh be used near the step comer for computational accuracy. It should be also noted that 
the law of the wall does not formally apply to separated turbulent boundary layers. However, since 
the separation point is fixed at the comer of the backstep — and the flowfield is solved iteratively 
so that the friction velocity w T can be updated until it converges — major errors do not appear to 
result from its use. 

First, results will be presented for the standard K-e model. For this case — as well as the other 
results to follow — computed results for the mean velocity streamlines, the streamwise mean ve- 
locity profiles, the streamwise turbulence intensity profiles and the turbulence shear stress profiles 
are compared with the Kim et a/. 30 experimental data as updated by Eaton and Johnston. 31 In fig- 
ure 2(a) the computed streamlines are shown, indicating reattachment at X^H =6.12 — a result 
which is approximately 15% underprediction of the experimental reattachment point of X^H = 7.1. 
In figure 2(b), the streamwise mean velocity profiles predicted by the standard K-e model are com- 
pared with the experimental data. Except in the vicinity of the reattachment point, the comparisons 
are fairly good. More serious discrepancies between the model predictions and the experimental 
data occur in the initial part of the recovery zone for the streamwise turbulence intensity profiles 
as shown in figure 3(a). However, the model predictions for the turbulence shear stress profiles are 
reasonably good as can be seen from figure 3(b). 

Next, we consider the computations based on the recursion-/?/VG model. For computational ef- 
ficiency, the nonlocal kernels g\ and g 2 are approximated by one-dimensional delta functions. This 
is appropriate if k c » 1 . In Fig. 4, the variation of k c (normalized by the step height and expressed 
as k c H) is shown at various locations downstream of the step. As can be seen, the magnitude of k c 
is at least an order of magnitude larger than the grid size employed (k c H > 1 for the 200 x 100 mesh 
used). Thus the local representation of g\ and g 2 is justified. The explicit functional form of the 
coefficients C R] and C R2 , (29), is much more difficult to evaluate. As a first attempt at applying 
this recurs ion -RNG model we make a lowest order approximation that these coefficients are con- 
stants. Clearly, especially in light of the great success in enforcing a functional ansatz on these 
types of coefficients, 9 this will underestimate the robustness and accuracy of the recursion-/? /VG 
formulation. Nevertheless, it is deemed more appropriate here to apply the theory without adding 
on these extraneous functional ansatzes. It is also evident that in applying the one-dimensional lo- 
cal form of g l (|x - x'|) and g 2 (|x - xj) , the angular dependences of the second term in (22) is 
much stronger than that of the first term in (22). This is because of the higher order derivatives 
present. These angular integrations are also dependent on spatial position due to the changing flow 
characteristics. Again, for simplicity, we approximate these angular integrations by constants, and 
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thus put the model to an even stricter test. In particular, for the present analysis, and are 
taken to be 0.025 and 0.0342 x 10 -3 , respectively. 

Now, it will be demonstrated that the use of the proposed recursive RAG-based model can yield 
a more significant improvement in the results. The computed streamlines for the flow field shown 
in figure 5(a) have a mean reattachment point, X r /H ~ 6.72, a result which is about 5% lower than 
the experimental value. The corresponding mean velocity profiles shown in figure 5(b) indicate 
very good agreement with the experimental results. It should be noted that computations performed 
based on the e -RNG - based model for the same flow conditions and model constants yield a mean 
reattachment point, X^H ~ 6.42 (not shown herein). The difference in the size of the separated flow 
region can be clearly attributed to the contributions from the Reynolds stress terms representing 
the local interactions effects in the recursion RNG model. 

Furthermore the overall agreement between the turbulence intensity and the shear stresses with 
the experimental data shown in figures 6(a) and 6(b) are also good. The most notable difference 
between the predictions of the RAG-based models and the standard K-e model lies in the stream- 
wise turbulence intensity (figure 6(a)). The slight trough shaped variation predicted in this region 
is consistent with more recent independent experiments. 35 

In addition, to illustrate the differences associated with the modeling of the Reynolds stresses, 
the variation of the turbulent eddy-viscosity normalized with respect to its molecular counterpart, 
v r / v 0 , is shown at several locations downstream of the backward-facing step. As can be seen, in 
the recirculation region the eddy viscosity predicted by the standard K-e model is generally larger 
than that due to the recursion RNG leading to substantial reduction in the size of the separated flow 
region (cf., figures 2(a) and 5(a)). 

The wall pressure coefficient is an important parameter for engineering applications. In figures 
8(a)-(b), the pressure coefficients C p (= 2 [p -p r ]/U r 2 , where p r and U r are the reference pressure 
and velocity which are taken at the centerline of the inlet) obtained from the standard and the con- 
ventional and the RNG - based K-e models at the top and bottom walls are compared with the exper- 
imental data of Eaton and Johnston. 31 As can be seen, both the standard and the RiVG-based K-e 
models perform comparably well in reproducing the experimental trends. The skin friction coeffi- 
cients Cy-= 2 u 2 lJJ r 2 obtained from the standard and recursion-RAG K-e models are compared 
with the scaled experimental data of Driver and Seegmiller 36 for the bottom wall in figure 8(c). 
Here, we make use of the fact that the ratio Cf/Cfa, when taken as a function of the normalized 
distance (X-X r )/X n is independent of the expansion ratio (given that is the fully-developed 
skin friction coefficient andX r is the reattachment point). As can be seen, the recursive RAG-based 
K-e model performs better, however, both models are probably within the uncertainty of the exper- 
imental data. 
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4. CONCLUSIONS 


The recursion renormalization group ( RNG ) theory is utilized to develop Reynolds stress clo- 
sure models for the prediction of turbulent separated flows. Unlike the small parameter € -RNG 
model that only consider the interaction of the long wavelength modes, the recursion-/? A/G model 
includes both the local and nonlocal interaction of all the relevant resolvable scales. A formal de- 
velopment of the Reynolds stress model is presented and the resulting higher order terms by con- 
sideration of the local interaction effects are shown to be qualitatively similar to those recently ob- 

78 

tained by considering a hierarchy of second-order closure models. 

The ability of the proposed model to accurately predict separated flows is analyzed from a com- 
bined theoretical and computational standpoint by considering turbulent flow past a backward fac- 
ing step as a test case. The results obtained based on detailed computations demonstrate that the 
proposed recursion-/? A/G model can yield very good predictions for the turbulent flow of an in- 
compressible viscous fluid over a backward- facing step. It should be remembered that the deficien- 
cies of two-equation models are well established, particularly in turbulent flows with body forces 
or Reynolds stress relaxation effects . 1 Consequently, the findings of this study should not be inter- 
preted as an unequivocal endorsement of two-equation RNG models. Nonetheless, this study 
shows that properly calibrated two-equation turbulence models, which account for the anisotropy 
of the turbulent stresses, can be effective for the prediction of turbulent separated flows. 
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(a) Streamlines 



(b) Dimensionless mean velocity profile 
( Computations; 

o Experiments of Kim etal, 1980; Eaton & Johnston, 1981) 


Figure 2 Computed mean flow field for the standard K-e model 
[E = 1 :3; Re = 132,000; 200x100 mesh] 
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Figure 3 Computed turbulence stresses for the standard K-e model [E = 1:3; 

Re = 132,000; 200x100 mesh; computations; o experiments of Kim et 

al, 1980; Eaton & Johnston, 1981] 
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Figure 4 Computed profiles of the cut-off wavenumber for the recursion 
RNG K-t model [E = 1:3; Re = 132,000; 200x100 mesh; — computations; 
o experiments of Kim et al, 1980; Eaton & Johnston, 1981] 
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(b) Dimensionless mean velocity profile 
( Computations; 

o Experiments of Kim etal, 1980; Eaton & Johnston, 1981) 


Figure 5 Computed mean flowfield for the recursion RNG K-t model 
[E = 1 :3; Re = 132,000; 200x100 mesh] 


- 24 - 




0.0 0.1 0.2 0.1 0.2 0.1 0.2 0.1 0.2 0.1 0.2 0.1 0.2 



(a) Turbulence intensity 



(b) Turbulence shear stress 


Figure 6 Computed turbulence stresses for the recursion RNG K-e model 
[E = 1:3; Re = 132,000; 200x100 mesh; — computations; o experiments of 
Kim et al, 1980; Eaton & Johnston, 1981] 
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Figure 7 Computed profiles of eddy-viscosity [E = 1:3; Re = 132,000; 
200x100 mesh; computations for the standard K-e model; com- 

putations for the recursion RNG fC-e model] 
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(b) Pressure coefficient along the top wall 



Figure 8 Comparison of the predicted wall distribution with experiments 
[ computations for the standard K-e model; computations for the re- 

cursion RNG K - 8 model; o experiments; E= 1:3; Re = 132,000] 
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